A First Principles Density-Functional Calculation of the Electronic and Vibrational 

Structure of the Key Melanin Monomers* 
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We report first principles density functional calculations for liydroquinone (HQ), indolequinone 
(IQ) and semiquinone (SQ). Tliese molecules are believed to be the basic building blocks of the 
eumelanins, a class of bio-macromolecules with important biological functions (including photo- 
protection) and with potential for certain bioengineering applications. We have used the ASCF 
(difference of self consistent fields) method to study the energy gap between the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), Ahl- We show 
that Ahl is similar in IQ and SQ but approximately twice as large in HQ. This may have important 
implications for our understanding of the observed broad band optical absorption of the eumelanins. 
The possibility of using this difference in Ahl to molecularly engineer the electronic properties of 
eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from 
first principles. Each of the molecules have significantly different infrared and Raman signatures, 
and so these spectra could be used in situ to non-destructively identify the monomeric content of 
macromolecules. It is hoped that this may be a helpful analytical tool in determining the structure 
of eumelanin macromolecules and hence in helping to determine the structure-property-function 
relationships that control the behaviour of the eumelanins. 

PACS numbers: 87.15.-v, 82.35.Cd, 87.64.Je 



I. INTRODUCTION 

The melanins are an important class of pigmen- 
tary macromolecule found throughout the biosphere 0. 
Pheomelanin (a cysteinyl-dopa derivative) and eumelanin 
(formed from 5-6-dihydroxyindolequinone and other in- 
dolequinones) are the predominant forms in humans, and 
act as the primary photoprotectant in our skin and eyes. 
Consistent with this role, all melanins show broad band 
monotonic absorption in the UV and visible in the range 
1.5 to 5 eV In contrast, other biomolecules such 
as proteins and nucleic acids show only well defined ab- 
sorption peaks around 280 nm (4.5 eV) and little ab- 
sorption below that They are also efficient free rad- 
ical scavengers and antioxidants In direct contradic- 
tion with these photoprotective properties, both pheome- 
lanin and eumelanin are irnplicated in the development 
of melanoma skin cancer |^. For this reason, the pho- 
tophysics, photochemistry and photobiology of melanins 
are subjects of intense scientific interest. 

Despite work over several decades (with respect to 
eumelanin in particular), the more general structure- 
property-function relationships that control the be- 
haviour of these important bio-macromolecules are still 
poorly understood 0. The melanins are difficult 
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molecules to study: they are chemically and photochem- 
ically stable, and arc virtually insoluble in most com- 
mon solvents. It is fairly well accepted that eume- 
lanins are macromolecules of the various redox forms 
of 5,6-dihydroxyindolequinone (DHI or HQ) and 5,6- 
dihydroxyindole 2-carboxylic acid (DHICA) P, . How- 
ever, major questions still remain concerning their ba- 
sic structural unit 01 ■ Two opposing schools of thought 
exist: i) that eumelanins are composed of highly cross- 
linked extended hetero-polymers based upon the Raper- 
Mason scheme and ii) that eumelanins are actually 
composed of much smaller oligomers condensed into 4 or 
5 oligomer nano-aggregates [g]- Clearly, this is a fun- 
damental issue, and is the starting point for the con- 
struction of consistent structure-property-function rela- 
tionships. 

The answer to this question also has profound implica- 
tions for our understanding of the condensedjphase prop- 
erties of melanins. In 1960 Longuet-Higgins f9l|, in a land- 
mark paper, proposed that many of the physical proper- 
ties of melanins could be understood if they were semi- 
conductors. This proposition was lent further support by 
Pullman and Pullman [lol | in 1964 who applied molecu- 
lar orbital theory in its simplest form (the Hiickel ap- 
proximation) to 5,6-dihydroxyindole and one particular 
dimer type. This theoretical work was followed in 1974 by 
the experimental observations of McGinness, Corry and 
Proctor who demonstrated that a pellet of melanin could 
be made to behave as an amorphous electrical switch 
PH . They postulated that these materials may indeed 
be disordered organic semiconductors consistent with the 
Longuet-Higgins theory, and the newly developed mod- 
els of amorphous inorganic semiconductors [T^l- Several 



studies since have also claimed to show that melanins 
in the condensed solid-state are semiconductors 0, Q| . 
However, it is by no means certain that the conductivity 
reported in any of these experimental studies is electronic 
in nature. A clear idea of the basic structural unit is crit- 
ical to developing a consistent model for condensed phase 
charge transport in such disordered organic systems. It is 
also important in the context of "molecularly engineer- 
ing" melanins to have the ability to create or enhance 
functionality in high technology applications such as bio- 
sensors and bio- mimetic photovoltaics [l5llT6l|. 

It is well known that disorder plays a crucial role in de- 
termining the charge transport properties of amorphous 
inorganic semiconductors. In addition to the usual im- 
purity related disorder, organic semiconductors (of which 
eumelanin may well be an exotic example), present sev- 
eral different sources of disorder associated with struc- 
tural heterogeneity. For example, if one considers ei- 
ther the hetero-polymer or oligomer models, there are 
many different ways of constructing a eumelanin macro- 
molecule based upon the monomer sequences, cross- 
linking positions or tautomer combinations. These pos- 
sibilities have been studied theoretically for DHI-based 
macromolecules using a number of quantum chemical 
techniques. Notably, Galvao and Caldas 0, 0, 
used the Hiickel approximation to construct model homo- 
polymers. Intriguingly, they found that many of the im- 
portant "semiconductor related properties of eumelanin" 
stabilise at a relatively small number of monomer units 
- five or six in fact. Although they did not consider 
inter-chain effects, and essentially treated chain termina- 
tions as end defects, their work was the first indication 
that large, extended hetero-polymeric structures are not 
required to explain the physical properties of melanins. 
Bolivar-Marinez et al. and Bochenek and Gudowska- 
Nowak Jl] have used the intermediate neglect of differ- 
ential overlap (INDO) and other semi-empirical methods 
to perform similar calculations for oligomers, whilst other 
authors have utilised the power of density functional the- 
ory to perform first principles calculations |2^ on sin- 
gle monomers or dimers. 

Il'ichev and Simon considered the hydroquinone, 
(HQ) and indolcquinone (IQ) forms (see figure^, and de- 
termined them to be the most stable and hence the most 
likely forms of the monomers. However, several pulse ra- 
diolysis studies [H H IH HI have indicated the 

presence of the semiquinone (SQ - also shown in figure 
nj, a tautomer of IQ, and so we also consider SQ here. 

Motivated by the remarkable optical absorption char- 
acteristics of eumelanins, many of these authors have fo- 
cused on simulating the absorption spectra of HQ, IQ, 
SQ and related small oligomers 0, |23, H^- However, 
knowledge of the electronic states around the highest 
occupied molecular orbital (HOMO) and lowest unoc- 
cupied molecular orbital (LUMO) levels is only part of 
the story. A great deal of information is also contained 
within the vibrational and rotational structure of the 
molecules. In particular, the IR and Raman spectra may 



be useful analytical tools. This has proven to be the 
case in other biomolecules 0. They can be used to con- 
firm the accuracy of the first principles predictions by 
direct comparison with experimental evidence, and also 
to provide useful insight into local environment phenom- 
ena such as solvent-solute interactions. Furthermore, our 
IR and Raman spectra calculations have the advantage 
of coming from first principles calculations and not the 
semi-empirical methods often employed to calculate the 
absorption spectra. 

In this paper we present a full density functional 
theory (DFT) analysis of the electronic, vibrational 
and rotational structure of the indolcquinone (IQ), the 
semiquinone (SQ) and the hydroquinone (HQ). Our 
study is a considerable extension of previously published 
work since we: i) present calculated IR and Raman spec- 
tra for the first time, and ii) use the difference of self con- 
sistent fields (ASCF) approach to gain a more accurate 
understanding of the LUMO level. We find that the elec- 
tronic properties of the molecules are highly dependent 
upon the redox form. Hence, we can say that oligomers 
or hetero-polymers consisting of one or more redox forms 
will have large variations in the local chemical potential, 
or, in the language of the tight binding model, strongly 
random site energies and possibly also random hopping 
integrals. This is to say that such hetero-structures will 
be highly disordered. 

Our work is motivated by the desire to understand 
the implications of the redox form of the basic structural 
unit of eumelanins on the materials bulk properties. Such 
knowledge is a key starting point in any attempts to un- 
derstand charge transport in and the optical properties 
of these disordered heterogeneous organic conductors. It 
also has profound implications in unravelling the myster- 
ies of melanin biological functionality, and attempts to 
molecularly engineer melanin-like molecules for techno- 
logical applications. These preliminary theoretical cal- 
culations are part of an ongoing quantum chemical, ex- 
perimental solid-state and spectroscopic program aimed 
at gaining a more complete understanding of melanin 
structure- property- function relationships. 



II. CALCULATION DETAILS 

The chemical and electronic structures and the Raman 
and IR spectra were found from first principles DFT cal- 
culations. We have performed our calculations using the 
Naval Research Laboratory Molecular Orbital Library 
(NRLMOL) ,2a^^ 32, 33, 34, 35). NRLMOL per- 
forms massively parallel electronic structure calculation 
using gaussian orbital methods. In particular, for a dis- 
cussion of the calculation of the Raman and IR spectra 
see Refs. [S^, 0, 113 • Throughout we have used the 
Perdew, Burke and Ernzerhof (PBE) exchange cor- 
relation functional, which is a generalised gradient ap- 
proximation (GGA) containing no parameters. For each 
molecule we have fully relaxed the geometry with no sym- 
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FIG. 1: Schematic representation of the 5,6- 
dihydroxyindolequinone (DHI) redox forms: hydroquinone 
(HQ) and the tautomers indolequinone (IQ) and semiquinone 
(SQ). The numbers correspond to those used in table llTI 



TABLE I: The total energy (in eV) of the two lowest energy 
tautomers of indolequinone (IQ); semiquinone (SQ) and 2- 
semiquinone (2-SQ) which is simply SQ with the hydrogen 
atom bonded to Oil rather than OlO (c.f. figure 0. The 
energies are quoted relative to the energy of IQ. In each case 
we fully relaxed the geometries as shown in figures |5] and 
^ Our results for HQ, IQ and SQ are consistent with the 
calculations of Il'ichev and Simon [2^ |. The concentrations at 
300 K are in direct contradiction to the experimental evidence 

m m m in in m that sq is one of the building blocks 

of eumelanin. 



Molecule 


Name in 
Ref. 22 


Energy (eV) 


Concentration 
at r = 300 K 


IQ 


3a 





87% 


2-SQ 


3d 


0.05 


13% 


SQ 


3c 


0.19 


< 0.1% 



TABLE II: 


The calculated bond lengths (in 


A) and bond 


angles (in de 


grees) are found to be 


in good af 


jreement with 


previous work [3, l20|. (Note that there is no 


experimental 


data to compare with.) The atom 


numbers 


correspond to 


those shown 


in figure 








IQ 


SQ 


HQ 


C1-C2 


1.423 


1.513 


1.586 


C2-C3 


1.388 


1.486 


1.458 


C3-C4 


1.408 


1.351 


1.366 


C4-C5 


1.421 


1.488 


1.483 


C5-C6 


1.401 


1.435 


1.361 


C6-C1 


1.392 


1.362 


1.457 


C4-C7 


1.435 


1.457 


1.443 


C7-C8 


1.376 


1.358 


1.365 


C8-N9 


1.385 


1.435 


1.390 


N9-C5 


1.382 


1.310 


1.381 


Cl-OlO 


1.375 


1.353 


1.226 


C2-011 


1.376 


1.226 


1.226 


C1-C2-C3 


120.69 


117.88 


118.32 


C2-C3-C4 


119.89 


119.61 


119.72 


C3-C4-C5 


118.45 


120.78 


121.03 


C4-C5-C6 


122.20 


121.32 


124.60 


C5-C6-C1 


118.18 


118.79 


118.34 


C6-C1-C2 


120.59 


121.62 


117.98 


C3-C4-C7 


134.47 


134.88 


132.73 


C6-C5-N9 


130.61 


127.05 


129.38 


C4-C7-C8 


107.08 


105.05 


107.23 


C7-C8-N9 


109.33 


113.75 


110.86 


C8-N9-C5 


109.32 


105.23 


109.65 


N9-C5-C4 


107.19 


111.63 


106.02 


C5-C4-C7 


107.08 


104.34 


106.24 


O10-C1-C6 


123.16 


125.59 


122.87 


O10-C1-C2 


116.24 


112.79 


119.15 


011-C2-C1 


116.15 


120.98 


118.89 


011-C2-C3 


123.16 


121.14 


122.79 



metry constraints. 

III. RESULTS AND DISCUSSION 

To benchmark our calculations we report the bond 
lengths and bond angles of HQ, IQ and SQ found in our 
calculations in table Our results are in good agree- 
ment with those calculated by Bolivar-Marinez et al. [20i | 
and Stark et al. j^. All three molecules are planar to 
within numerical accuracy (c.f., [EE^)- 

We find, in agreement with Il'ichev and Simon [2^ . 
that SQ is not the lowest energy tautomer of IQ (see table 
P) . Indeed we find that at 300 K there should be less than 
0.1% SQ. Although we should stress that this calculation 
is for a single molecule in vacuo it is, of course, in direct 
contradiction to the evidence jH IH IH IH lllllll that 
SQ is one of the building blocks of eumelanin. However, 
since SQ is widely thought to play an important role in 
eumelanin, we also consider it here. 

As well as the ground state geometries shown in figures 



|21E|and^we also found stable geometries of HQ and SQ 
with the H-O-C bond angle increased by almost 180°. 
These alternative structures have significantly higher en- 
ergies than the ground states. However, it is possible that 
the high energy structures may be stabilised by hydrogen 
bonding in a polar solvent. We therefore also considered 
the system SQ-I-6H2O. We found that in this system the 
H-O-C bond is slightly larger than in SQ in vacuo (see 
figure E)). We also found a significant hydrogen bonded 
network around the H-0-C-C=0 group which is dragged 
out of the plane of the molecule by interactions with the 
water molecules. This hydrogen bonding network is prob- 
ably related to the fact that eumelanins are only soluble 
in polar solvents. 

The Hiickel model calculations of Galvao and Caldas 
ITtI indicate that the HOMO-LUMO gap (Ahl) niay be 
significantly larger in HQ than in either IQ or SQ. This 
has important implications in the context of understand- 
ing the basic eumelanin structural unit. Fundamentally, 
a combination of HQ and either IQ or SQ tautomers 
are present in an ensemble. Whether that ensemble be 
oligomeric or polymeric it can be expected to lead to an 
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apparent broadening of the absorption profile. In this 
situation, the key question is: how many monomers of 
HQ, IQ and SQ are required to create the broad band 
UV and visible absorbance over the range 1.5 to 5 eV 
that is macroscopically observable for eumelanins? The 
answer to this question may well allow us to assess the 
feasibility of the various structural models. Furthermore, 
when one considers the possibility that molecularly en- 
gineered forms of eumelanin may be useful as functional 
materials in electronic devices and sensors 0, , a crit- 
ical design parameter is likely to be the "semiconductor" 
gap. It has previously been shown that, at least 
within the Hiickel approximation for homopolymers, the 
HOMO-LUMO gap of the monomer is closely related to 
the HOMO-LUMO gap of the infinite polymer, which, 
for the polymer model in the condensed solid state, cor- 
responds to the semiconducting gap provided there is 
only weak coupling between polymers. Thus controlling 
the HOMO-LUMO gap of macromolecules may provide 
a route to controlling the semiconducting gap of eume- 
lanins. This suggestion clearly needs further investiga- 
tion, but controlling the HQ to IQ or SQ ratio may be 
a possible route to this form of property manipulation. 
An important first step to investigating either of these 
proposals is the calculation of the HOMO-LUMO gaps 
of the molecules by more reliable methods than Hiickel 
theory. 

In table Uni we compare the HOMO-LUMO gap found 
in our calculations with that found from the Hiickel 
method For comparison with the Hiickel calcula- 

tions we need to estimate the hopping integral, /?, in the 
Hiickel Hamiltonian, 

ii = + a'5^^3)c\„c^„ - "^Lcj-t, (1) 

where a is the on-site energy, a' is the difference in the 
on-site energy between a nitrogen and a carbon atom (wc 
take the labelling convention from figure ^ ensuring that 

the only nitrogen is in position i = 9), c^P annihilates 
(creates) an electron with spin a on site i and [ij) indi- 
cates that the sum is over nearest neighbours only. For 
hydrocarbons j3 is typically of order 2.5 eV 37]. Using 
this value of (3 we find that the HOMO-LUMO gaps for 
HQ, IQ and SQ found from the Hiickel approximation 
are consistent with those found in our DFT calculations 
using the FEE functional. However, because DFT is a 
theory of the ground state, these calculations represent 
the energy gap between Kohn-Sham eigenvalues and not 
the true HOMO-LUMO gap of the molecules. This is 
known as the band gap problem |38j| . Additionally, it is 
accepted that the PEE functional can significantly under- 
estimate the HOMO-LUMO gap. Therefore we have also 
employed the ASCF method [33 to calculate Ahl- Our 
results clearly reproduce the trends seen in the time de- 
pendent density functional theory (TDDFT) calculations 
of Il'ichev and Simon |23| for IQ and SQ. The difference 
between the ASCF and TDDFT results for HQ probably 
indicates that the true value of Ahl is intermediate to 



TABLE III: The HOMO-LUMO (highest occupied molecular 
orbital-lowest unoccupied molecular orbital) gap, Ahl, in eV. 
The HOMO-LUMO gap found from both the ASCF method 
and from a simple density functional calculation (both cal- 
culations use the PBE functional [s^) are compared with 
that found in the Hiickel approximation (after Il7l|, t aking 
(3 = 2.5 eV which is the correct order of magnitude |37ll . c.f., 
Q). It can be seen that the HOMO-LUMO gap of both the 
indolequinone (IQ) and the semiquinone (SQ) is significantly 
underestimated by both the Hiickel approximation and the 
simple interpretation of the Kohn-Sham eigenvalues. How- 
ever, these methods find a HOMO-LUMO gap for the hydro- 
quinone (HQ) that is more consistent with the, more reliable, 
ASCF method. The results of the ASCF method indicate 
that the HOMO-LUMO gap is approximately the same in 
IQ and SQ but about a factor of two larger in HQ. These 
results correctly reproduce the trends seen in the more accu- 
rate, more more computationally expensive time dependent 
density functional theory (TDDFT) calculations of Il'ichev 
and Simon |23 |. This is consistent with the threshold for op- 
tical absorption found from semi-empirical methods 5, 20]. 
The large difference in the HOMO-LUMO gap between HQ 
and the other molecules may be important for "molecular 
engineering" applications and for explaining the broad band 
optical absorption of eumelanin. In particular, if the model 
of the eumelanins as amorphous semiconductors is correct it 
may allow the chemical tuning of the semiconducting band 
gap (c.f. Ref. 18]). 





IQ 


SQ 


HQ 


TDDFT (PBE/B3LYP) [22] 


1.82 


1.50 


4.53 


TDDFT (B3LYP) |22] 


1.79 


1.43 


4.30 


NRLMOL (ASCF/PBE) 


2.02 


1.12 


3.61 


NRLMOL (PBE) 


1.07 


0.80 


3.48 


Hiickel [17] 


1.3 


0.84 


3.4 



the values predicted by the two different methods as it 
is extremely unlikely that either method is incorrect by 
a large enough margin to allow the other to be correct. 

The calculated ASCF HOMO-LUMO gap is also con- 
sistent with previously published semiempirical ZINDO 
(Zerner's intermediate neglect of differential overlap) cal- 
culations of the optical absorption spectra |^ |23]. In 
particular Eolivar-Marinez et al. noted that "in the neu- 
tral state the threshold for optical absorption is around 
2.0 eV for the IQ and SQ, while in the case of the HQ it is 
roughly 3.8 eV." Taking these numbers as semiempirical 
estimates of the HOMO-LUMO gap they are in excellent 
agreement with our ASCF calculations (c.f. table IIIl|l . 

It is also interesting to note that the change in the SQ 
conformation observed in the calculations for the SQ -I- 
6H2O system caused a 7% decrease in Ahl relative to 
SQ in vacuo. It is possible that in a polar solvent a range 
of conformations exist which may lead to a broadening 
of the optical absorption ^ . 

In terms of Hiickel theory, the difference between Ahl 
for HQ and IQ/SQ can be understood in terms of the 
level of delocahsation of the HOMOs and LUMOs. It 
should be noted that the aromatic ring in HQ contains 
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FIG. 2: The electron density in the highest occupied molecu- 
lar orbital (HOMO) (top) and the lowest unoccupied molecu- 
lar orbital LUMO (bottom) of hydroquinone (HQ). The atoms 
are colour coded as follows: carbon - green, nitrogen - blue, 
oxygen - red and hydrogen - white. The LUMO electron den- 
sity is in good agreement with the semiempirical calculations 
of Boliar-Marinez et al. I2Q|, however there is some discrep- 
ancy in the HOMO electron density. 



a resonating valence bond, whereas this feature is absent 
in IQ and SQ. This can be contrasted with the Hiickel 
model result for benzene which predicts Ahl — 2/3 for 
the neutral molecule but if one double bond in the ring is 
replaced by a single bond, as in, for example, the reduced 
form CgHg, then Ahl ~ 0. On this basis one might well 
expect the HOMO-LUMO gap of HQ to be significantly 
larger than those of IQ and SQ. However, detailed plots 
of the electron densities in the HOMO and LUMO of 
the three molecules (figures |21 El and Q indicate that 
the real situation is somewhat more complicated. The 
electron densities are significantly more localised in the 
HOMOs of IQ and SQ than the electron density of the 
HOMO of HQ is. In particular, the electron densities in 
the HOMOs and LUMOs of IQ and SQ are not as similar 
as one might expect from figure ^ However, even in the 
ASCF calculation, the HOMO-LUMO gap is a factor of 
two larger in HQ than it is in SQ and IQ. 

In the absence of knowledge of the chemical structure 
of the eumelanins it is difficult to predict how a change 
in the HOMO-LUMO gap of the monomer will effect the 
electronic structure. However it is likely that changing 
Ahl by a factor of two will have a dramatic effect on 
the electronic structure of the melanin macromoleculc. 
For example, it is reasonable to expect that the HOMO- 
LUMO gap of a macromolecule will be related to Ahl 
(c.f. Ref. 18]). Thus we expect that the HOMO-LUMO 
gap for a macromolecule of HQ will be significantly larger 
than that of a macromolecule of IQ or SQ. For a macro- 




FIG. 3: The electron density in the highest occupied molecu- 
lar orbital (HOMO) (top) and the lowest unoccupied molecu- 
lar orbital LUMO (bottom) of indolequinone (IQ). The atoms 
are colour coded as follows: carbon - green, nitrogen - blue, 
oxygen - red and hydrogen - white. The LUMO electron den- 
sity is in good agreement with the semiempirical calculations 
of Boliar-Marinez et al. (20( . however there is some discrep- 
ancy in the HOMO electron density. 




FIG. 4: The electron density in the highest occupied molecu- 
lar orbital (HOMO) (top) and the lowest unoccupied molecu- 
lar orbital LUMO (bottom) of semiquinone (SQ). The atoms 
are colour coded as follows: carbon - green, nitrogen - blue, 
oxygen - red and hydrogen - white. Both the HOMO and 
LUMO electron densities are in good agreement with the 
semiempirical calculations of Boliar-Marinez et al. |2(i^ . 
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TABLE IV: The predicted Phonon (IR) and Raman spectra 
for hydroquinone (HQ). 



Isotropic Raman Total Raman 
Frequency IR intensity scattering scattering 
(cm~^) (D^/amuA^) activity activity 

(AVamu) (AVamu) 




FIG. 5: The stable structure of semiquinone (SQ) and 6 II2O 
molecules. Note the hydrogen bonded networks around the 
H-0-C-C=0 group, which is dragged slightly out of the plan 
by the interactions with the solvent. The H-O-C bond angle 
has increased relative to that of the SQ molecule in vacuo. 



molecule containing both HQ and, say, IQ this change in 
the electronic structure will act as a source of disorder 
and thus have a dramatic effect on the electronic trans- 
port properties. 

It is also difficult, without detailed knowledge of the 
chemical structure of the eumelanins, to predict exactly 
how the variation of Ahl affects the optical absorption. 
However, it seems possible that the range of Ahl be- 
tween different molecules may produce the continuum of 
HOMO-LUMO gaps in eumelanin macromolecules and 
thus play an important role in explaining the observed 
the broad band optical absorption. 

It is therefore important to be able to identify the 
monomeric content of a sample of eumelanin, ideally this 
should be done in situ and non-destructively. Motivated 
by this fact we have calculated both the Raman and in- 
frared (IR) spectra of IQ, SQ and HQ, (see tables HVl Ivl 
and|Vl|)- This has several uses. Firstly, the prediction 
allows for experimental testing of the accuracy of our 
calculations. Secondly, it can be seen that there are no- 
table differences in both the Raman and IR spectra of the 
three monomers (see, for example, figures El and I?)) . The 
Raman and IR spectra could therefore be used for the 
in situ, non-destructive identification of the monomeric 
content of macromolecules. This could be particularly 
valuable in the engineering of devices from eumelanins, 
as, by varying the ratio of indole-quinones, it may be 
possible to achieve control of the band-gap of the ma- 
terial. This may also be a useful analytical tool for the 
analysis of the structure of both synthetic and naturally 
occurring melanin, and thus be helpful in determining 
the structure-property-function relationships that con- 
trol their behaviour. 

Comparison of our calculated IR and Raman spectra 
with standard tables for organic molecules ,41j shows 
that the calculated phonon spectra are in broad agree- 
ment with known IR absorption band values (see table 
IVII|) ■ It can also be seen from table IVlTl that the dif- 
ferences in the spectra, specifically which vibrations are 
absent for a given molecule can be understood in terms 
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TABLE V: The predicted Phonon (IR) and Raman spectra 
for indolcquinonc (IQ). 
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TABLE VI: The predicted Phonon (IR) and Raman spectra 
for semiquinone (SQ). 
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of their structure, c.f. figure^ as expected. 

It is important to stress that our calculated IR and 
Raman spectra are for gaseous phase monomers. Thus, 
the comparison of such results with experiment may be 
complicated by several factors. Firstly the formation of 
macromolecules may result in a significant change in the 
spectra. However, it is known that in other biomolecules 
the effect of polymerisation only slightly broadens the 
IR and Raman spectra, for example the spectra of sim- 
ple amino acids is only slightly broadened in proteins 
p^. Secondly solvent effects may largely 'wash out' the 



individual features of the spectra [s^ 0]. It is clear 
from our study of the SQ + 6H2O system that there will 
be significant hydrogen bonding in the solvent, hydrogen 
bonding can shift the phonon frequencies Therefore 
it remains to be seen if, experimentally, the resolution is 
good enough to resolve the differences in the spectra of 
the monomers. However, an obvious way to avoid solvent 
effects is to conduct the experiments in the solid state. 
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TABLE VII: Broad assignments of the strongest bands in the simulated phonon (IR) spectra of HQ, SQ and IQ (FIG. 6). The 
band assignments were made according to standard tables |40l. |41|. 



Band (cm ^) 



Group /Vibration 



HQ 



SQ 



IQ 



-1100 

-1340 

-1450-1590 

-1620-1650 

-3590-3600 

-3700 



V (C-OH) strong 

V (Aromatic C-N) strong 

V (Aromatic C-C) strong 

V (quinone C=0) strong 

V (N-H) medium 
V (O-H) strong 



present 
present 
present 
absent 
present 
present 



present 
absent 
present 
present 
absent 
present 



absent 
present 
present 
present 
present 
absent 



IV. CONCLUSIONS 

We have carried out first principles density func- 
tional calculations for the hydroquinone (HQ), the in- 
dolequinone (IQ) and the semiquinone (SQ). The calcu- 
lated gaseous phase structure is in good agreement with 
previous calculations. We have used the ASCF method 
to study the HOMO-LUMO gap. Our results are con- 
sistent with the threshold for optical absorption found 
by semiempirical methods. Specifically, we found that 
the HOMO-LUMO gap is similar in IQ and SQ but ap- 
proximately twice as large in HQ. The possibility of us- 
ing this difference in the HOMO-LUMO gap to engineer 
the electronic properties of eumelanins at the molecular 
level has been discussed. We have also suggested that 
the difference in the HOMO-LUMO gap of the different 
monomers could lead to a large range of HOMO-LUMO 
gaps in eumelanin macromolecules and thus be related 
to the observed broad band optical absorption. 

As the structure of macromolecules of these monomers 
is not known we have also calculated the IR and Raman 
spectra of the three monomers from first principles. A 
comparison of these results with experiment would repre- 
sent a stringent test of the density functional calculations 
as these calculations do not rely on additional semiem- 
pirical information. It was shown that the IR and Raman 



spectra have potential as a analytical tool for these mate- 
rials. Each of the monomers have significantly different 
spectral signatures, therefore the IR or Raman spectra 
could be used in situ to non-destructively investigate the 
monomeric content of macromolecules. It is hoped that 
this may be helpful in determining the structure of both 
natural and synthetic eumelanins, and hence aid our at- 
tempts to understand their biological functionality. 



Acknowledgements 

We would like to thank John Dobson, Joel Gilmore, Pe- 
ter Innis, Urban Lundin, Jeffrey Reimers, Jennifer Riesz 
Cathy Stampff and Linh Tran for helpful discussions. 
The work at the University of Queensland was supported 
by the Australian Research Council. BJP was supported 
in part by US Navy Grant N00014-03-1-4115. KB was 
supported in part by US Navy Grant N00014-03-1-4116. 
TB was supported in part by US Navy Grant N00014- 
02-1-1046. MRP was supported in part by ONR and 
the DoD HPC CHSSI initiative. BJP, RHM and MRP 
appreciate the hospitality of the Kavli Institute of Theo- 
retical Physics, where this work was initiated. The KITP 
is supported by NSF through grant PHY99-07949. 



[1] G. Prota, Melamns and Melanogenesis (Academic Press, 

New York, 1992). 
[2] M. L. Wolbarsht, A. W. Walsh, and G. George, Appl. 

Opt. 20, 2184 (1981). 
[3] K. E. van Holde, W. C. Johnson, and P. S. Ho, Princi- 
ples of Physical Biochemistry (Prentice-Hall, Englewood 

Cliffs, NJ, 1998). 
[4] H. Z. Hill, Melanin: Its Role m Human Photoprotection 

(Valdenmar Press, Overland Park, KS, 1995). 
[5] K. B. Stark, J. M. Gallas, G. W. Zajac, and J. T. Golab, 

J. Phys. Chem. B 107, 3061 (2003). 
[6] S. Ito, Biophys. Acta 883, 155 (1986). 
[7] Z. W. Zajac, J. M. Gallas, J. Cheng, M. Eisner, S. C. 

Moss, and A. E. Alvarado-Swaisgood, Biochim. Biophys. 

Acta 1199, 271 (1994). 
[8] C. M. R. Clancy, J. B. Nofsinger, R. K. Hanks, and J. D. 

Simon, J. Chem. Phys. B 104, 7871 (2000). 
[9] H. C. Longuet-Higgins, Arch. Biochem. Biophys. 88, 231 



(1960). 

[10] A. Pullman and B. Pullman, Biochim. Biophys. Acta 54, 
384 (1961). 

[11] J. McGinness, J. P. Corry, and P. Proctor, Science 183, 
853 (1974). 

[12] N. F. Mott and E. A. Davis, Electronic Processes in Non- 
Crystalline Materials (Clarendon Press, Oxford, 1971). 

[13] P. R. Crippa, V. Cristofoletti, and N. Romeo, Biochim. 
Biophys. Acta 538, 164 (1978). 

[14] M. M. Jastrzebska, V. Cristofoletti, and N. Romeo, 
Biochim. Biophys. Acta 7, 781 (1996). 

[15] P. Meredith and J. Riesz, cond- mat/0312277 (to appear 
in Photochem. Photobiol.). 

[16] J. Riesz, Honours thesis, Unversity of Queensland (2003). 

[17] D. S. Galvao and M. J. Caldas, j'. Chem. Phys. 88, 4088 
(1988). 

[18] D. S. Galvao and M. J. Caldas, J. Chem. Phys. 92, 2630 
(1990). 



8 



[19] D. S. Galvao and M. J. Caldas, J. Chem. Phys. 93, 2848 

(1990). 

[20] L. E. Bolivar-Marincz, D. S. Galvao, and M. J. Caldas, 

J. Phys. Chem. B 103, 2993 (1999). 
[21] K. Bochenek and E. Gudowska-Nowak, Chem. Phys. 

Lett. 373, 523 (2003). 
[22] Y. V. Il'ichev and J. D. Simon, J. Phys. Chem. B 107, 

7162 (2003). 

[23] C. Lambert, J. N. Chacon, M. R. Chedekel, E. J. Land, 

P. A. Riley, A. Thonpson, and T. G. Truscott, Biochem. 

Biophys. Acta 993, 12 (1989). 
[24] C. Lambert, E. J. Land, P. A. Riley, and T. G. Truscott, 

Biochem. Biophys. Acta 1035, 319 (1990). 
[25] A. T. Al-Kazwini, P. O'Neill, G. E. Adams, R. B. Cun- 

dall, B. Jacquet, G. Lang, and A. Junino, J. Phys. Chem. 

94, 6666 (1990). 
[26] A. T. Al-Kazwini, P. O'Neill, G. E. Adams, R. B. Cun- 

dall, G. Lang, and A. Junino, J. Chem. Soc, Perkin 

Trans. 2 12, 1941 (1991). 
[27] A. T. Al-Kazwini, P. O'Neill, R. B. Cundall, G. E. 

Adams, A. Junino, and J. Maignan, Tetrahedron Lett. 

33, 3045 (1992). 
[28] A. T. Al-Kazwini, P. O'Neill, G. E. Adams, R. B. Cun- 
dall, A. Junino, and J. Maignan, J. Chem. Soc, Perkin 

Trans. 2 13, 657 (1992). 
[29] M. R. Poderson and K. A. Jackson, Phys. Rev. B 41, 

7453 (1990). 

[30] K. A. Jackson and M. R. Pederson, Phys. Rev. B 42, 



3276 (1990). 

[31] M. R. Pederson and K. A. Jackson, Phys. Rev. B 43, 
7312 (1991). 

[32] A. A. Quong, M. R. Pederson, and J. L. Feldman, Solid 

State Commun. 87, 535 (1993). 
[33] D. V. Porezag and M. R. Pederson, Phys. Rev. B 54, 

7830 (1996). 

[34] D. V. Porezag, Ph.D. thesis, Technische Universitat 
(1997), http://archiv.tu-chemnitz.do/pub/1997/0025. 

[35] A. Briley, M. R. Pederson, K. A. Jackson, D. C. Patton, 
and D. V. Porezag, Phys. Rev. B 58, 1786 (1998). 

[36] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[37] J. P. Lowe, Quantum Chemistry (Academic Press, Lon- 
don, 1978). 

[38] R. O. Jones and O. Gunnarsson, Rev. Mod. Phys. 61, 
689 (1989). 

[39] A. B. Myers, Ann. Rev. Phys. Chem. 49, 267 (1998). 

[40] K. Nakanishi and P. H. Solomon, Infrared Absorption 
Spectroscopy (Holden-Day Inc., San Francisco, 1977). 

[41] E. B. Schrader, Infrared and Raman Spectroscopy, Meth- 
ods and Applications (VCH Publishers, New York, 1995). 

[42] A. Xie, L. van der Meer, W. Hoff, and R. H. Austin, 
Phys. Rev. Lett. 84, 5345 (2000). 

[43] C. Reichaxdt, Solvents and Solvent Effects in Organic 
Chemistry (VCH Verlagsgesellschaft, Basel, 1998). 



9 




FIG. 6: Detail of the phonon (IR) spectra of hydroquinone (HQ) (solid line), semiquinone (SQ) (dashed line) and indolequinone 
(IQ) (dotted Une). Clear differences can be seen between the spectra of the three monomers. For example, there are strong 
features around llOOcm"^ in both HQ and SQ but not in IQ. On the other hand, there are strong bands between ~1620 and 
1650cm~^ in SQ and IQ which are absent in HQ (for broad assignments see table IVm . This confirms that the IR spectra 
could be useful as in situ, non-destructive probes for identifying the monomeric content of eumelanins. To allow us to plot the 
spectra and to aid comparison with experiment we have broadened each peak by 6cm~^, a value that is comparable to that 
seen experimentally for other biomolecules in solution _3|. 




V (1/cm) 



FIG. 7: Detail of the Raman spectra of hydroquinone (HQ) (solid line), semiquinone (SQ) (dashed line) and indolequinone 
(IQ) (dotted line). Once again, there are clear differences between the spectra of the three monomers confirming that Raman 
spectra could be useful as in situ, non-destructive probes for identifying the monomeric content of eumelanins. To allow us to 
plot the spectra and to aid comparison with experiment we have broadened each peak by 6cm~^, a value that is comparable 
to that seen experimentally for other biomolecules in solution |^. 
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